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AUTOMATED COMPUTER-IMPLEMENTED METHOD 
AND SYSTEM FOR REORIENTING 
EMISSION COMPUTER TOMOGRAPHIC 
MYOCARDIAL PERFUSION IMAGES 

5 BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

One aspect of the present invention relates to an automated computer- 
implemented method and system reorienting emission computed tomographic (ECT), 
for example, single photon emission computer tomography (ECT) or position 
10 emission tomography (PET) heart images. More specifically, this invention relates 
to a method and system for aligning and centering ECT myocardial perfusion 
images. 

2. Background Art 

Emission computed tomographic (ECT) cardiac imaging can be used 
15 for the diagnosis of patients with coronary heart diseases (CHD), for example, 
patients with myocardial infarction or angina pectoris. An ECT heart study can be 
utilized to determine whether a patient's heart muscle or myocardium is healthy, 
diseased or scarred. In a typical study, two ECT image sets of a patient are 
acquired; the first after injection of a radiotracer of heart blood flow while the 
20 patient is at rest; the second after injection while the patient is performing a 
treadmill exercise stress test. If a particular region of the heart demonstrates 
reduced blood flow in both the rest and stress image sets, then it is likely that the 
heart muscle in that region is secured. If a particular region of the heart displays 
normal blood flow at rest but reduced blood flow in the stress image, this indicates 
25 that the coronary artery providing blood Co that is at least partially blocked, 
narrowed or stenosed. If a particular region of the heart is bright or normally 
perfused in both the rest and stress images, it is likely that the blood supply is 
normal and that particular region or the heart is healthy. 
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Typically, two-dimensional ECT images (tomographic slices) are 
generated orthogonal to (short axis slices) and parallel to (horizontal and vertical 
long axis slices) the heart's left ventricular (LV) long axis. Figures la-Id depict a 
three-dimensional rendering of the left ventricle and the typical slicing of the left 
5 ventricle into short axis and horizontal end vertical long axis sections. Figure la 
depicts a three-dimensional rendering of the LV within a coordinate system. A short 
axis slice is shown in Figure lb, a horizontal slice is shown in Figure lc, and a 
vertical long axis slice is shown in Figure Id. While a short axis slice, as depicted 
in Figure lb, is generally circular, the long axis slices, as depicted in Figures lc and 

10 Id are generally horseshoe shaped. Short axis two-dimensional images can be 
generated at discrete positions along the LV long axis between the base and apex of 
the heart, as depicted in Figure la. Horizontal long axis images can be generated 
at discrete positions along the short axis of the heart between the anterior and 
posterior walls of the heart. Vertical long axis images can be generated at discrete 

15 positions along the short axis of the heart between the septal and lateral heart walls. 

The optimal diagnostic accuracy of ECT myocardial perfusion studies 
requires the proper orientation, i.e. exact alignment, of the ECT images for visual 
and quantitative comparisons of one study to the other. Misalignment of the walls 
of the left ventricle between serial studies (e.g. stress and rest), can cause false 

20 diagnoses. Misalignment is commonly caused by failure to exactly position the 
patient in the ECT scanner for each of the two or more studies. In large and/or 
acutely ill patients this can be virtually impossible to accomplish. Misalignment also 
commonly occurs during the processing of the studies into tomographic slices. 
There have been many prior art proposals for minimizing the image artifacts caused 

25 by an inaccurate registration and alignment of patient studies. 

The approaches proposed thus far usually necessitate manual and 
interactive efforts to address misalignment. In light of the foregoing, an automated 
method and system for realigning and centering ECT myocardial perfusion images 
is needed. What is needed is a computer-implemented method or system for 
30 automatic alignment and centering of ECT cardiac images about the heart's left 
ventricular (LV) long-axis. 
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SUMMARY OF THE INVENTION 

One aspect of the present invention relates to a computer-implemented 
method and system for realigning and centering of ECT myocardial perfusion 
images. The present invention also includes a computer-implemented method and 
5 system for automatic realignment and centering of SPECT myocardial perfusion 
images about the left ventricle long-axis of the heart. One object of the present 
invention is to enhance reproducibility of quantitative measured ventricular 
perfusion and function. 

According to one embodiment of the present invention, an automated 

10 computer-implemented method for reorienting ECT myocardial perfusion images of 
a left ventricle (LV) of a heart is disclosed. The method includes receiving 
variously oriented tomographic images of a reconstructed volume of a heart; 
determining or receiving the LV long-axis and LV center (x c ,y c ,z c ) and determining 
or receiving axial limits (z apex and z base ) of the LV based on the tomographic images; 

15 determining or receiving the endocardial surface of the LV based on the 
tomographic images; determining a reorientation slice range based on the center and 
axial limits of the LV; receiving or reconstructing a number of slices (N) within the 
reorientation slice range; for each slice i=l to N, determining a center coordinate 
x[i], y[i] based on the endocardial surface and the area of the slice i within a 

20 reorientation slice range coordinate system; determining translation values Ax t , A^ 
and rotation values 0 X , 0 y based on the center coordinates x[i= 1 to N], y[i= 1 to N] 
to reorient the long axis of the LV to the z-axis and its origin of a reference 
Cartesian coordinate system; and automatically reorienting and realigning the 
tomographic images based on the translation values and rotation values to obtain 

25 reoriented ECT images reoriented and realigned to the long axis of the left ventricle. 

The reorienting step can include realigning, co-registering and 
centering the variously oriented input tomographic images to obtain a plurality of 
oriented ECT images. The reorienting step can be performed using an affine 
transform. The method can further include determining functional parameters based 
30 on the reoriented ECT images. The method can further include determining 
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perfusion parameters based on the reoriented ECT images. The variously oriented 
tomographic images are generated by perfusion imaging software. 

According to another embodiment of the present invention, a 
computer-implemented system for reorienting ECT myocardial perfusion images of 
5 a left ventricle (LV) of a heart is disclosed. The system includes one or more 
computers, configured to: receive variously oriented tomographic images of a 
reconstructed volume of a heart; determine or receiving the LV long-axis and LV 
center (x^y^zj and determine or receive axial limits (z apex and z base ) of the LV based 
on the tomographic images; determine or receive the endocardial surface of the LV 

10 based on the tomographic images; determine a reorientation slice range based on the 
center and axial limits of the LV; receive or reconstruct a number of slices (N) 
within the reorientation slice range; for each slice i=l to N, determine a center 
coordinate x[i], y[i] based on the endocardial surface and the area of the slice i 
within a reorientation slice range coordinate system; determine translation values 

15 AXj, Ayj and rotation values 0 X , 0 y based on the center coordinates x[i=l to N], 
y[i= 1 to N] to reorient the long axis of the LV to the z-axis and its origin of a 
reference Cartesian coordinate system; and automatically reorient and realign the 
tomographic images based on the translation values and rotation values to obtain 
reoriented ECT images reoriented and realigned to the long axis of the left ventricle. 

20 The one or more computers can be configured to realign, co-register 

and center the variously oriented input tomographic images to obtain a plurality of 
oriented ECT images. The one or more computers can be further configured to 
reorient and realign the tomographic images using an affme transform. The one or 
more computers can be further configured to determine functional parameters based 

25 on the reoriented ECT images. The one or more computers can be further 
configured to determine perfusion parameters based on the reoriented ECT images. 
The variously oriented tomographic images are generated by perfusion imaging 
software. 

According to yet another embodiment of the present invention, a 
30 computer-implemented apparatus for reorienting ECT myocardial perfusion images 
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of a left ventricle (LV) of a heart is disclosed. The apparatus includes means for 
receiving variously oriented tomographic images of a reconstructed volume of a 
heart; means for determining or means for receiving the LV long-axis and LV center 
(x c ,y c ,z c ) and means for determining or means for receiving axial limits (z apex and 
5 z base ) of the LV based on the tomographic images; means for determining or means 
for receiving the endocardial surface of the LV based on the tomographic images; 
means for determining a reorientation slice range based on the center and axial limits 
of the LV; means for receiving or means for reconstructing a number of slices (N) 
within the reorientation slice range; for each slice i = 1 to N, means for determining 

10 a center coordinate x[i], y[i] based on the endocardial surface and the area of the 
slice i within a reorientation slice range coordinate system; means for determining 
translation values AXj, and rotation values 0 X , 0 y based on the center coordinates 
x[i= 1 to N], y[i= 1 to N] to reorient the long axis of the LV to the z-axis and its 
origin of a reference Cartesian coordinate system; and means for automatically 

15 reorienting and realigning the tomographic images based on the translation values 
and rotation values to obtain reoriented ECT images reoriented and realigned to the 
long axis of the left ventricle. 

The means for reorienting can include a means for realigning, co- 
registering and centering the variously oriented input tomographic images to obtain 

20 a plurality of oriented ECT images. The means for reorienting can include a means 
for performing an affine transform. The apparatus can further include a means for 
determining functional parameters based on the reoriented ECT images. The 
apparatus can further include a means for determining perfusion parameters based 
on the reoriented ECT images. The variously oriented tomographic images are 

25 generated by perfusion imaging software. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The features of the present invention which are believed to be novel 
are set forth with particularity in the appended claims. The present invention, both 
30 as to its organization and manner of operation, together with further objects and 
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advantages thereof, may best be understood with reference to the following 
description, taken in connection with the accompanying drawings which: 

Figures la-Id depicts three-dimensional renderings and slices of a 
patient's left ventricle (LV); 

Figures 2a-2c depicts orthogonal slice views, i.e. ventricle long axis, 
mid short axis, and horizontal long axis, respectively, of the LV. The LV center 
of this particular example is X:33.0; Y:33.0 and Z:13.0; 

Figure 3 depicts a flowchart of a reorienting process according to one 
embodiment of the present invention; 

Figure 4 depicts a splash display of images (short axis (SA), 
horizontal long axis (HLA), and vertical long axis (VLA)) with various rotational 
misalignments; 

Figure 5 depicts a splash display of images (SA, HLA, and VLA) that 
have been auto-realigned about the long axis of the heart; 

Figure 6 depicts a graph showing the mean change in the functional 
estimates when the auto-realignment is employed; 

Figure 7 shows a splash display of images (SA, HLA, and VLA) with 
various misalignments and corresponding defect blackout maps and regional defect 
extent values; 

Figure 8 shows a splash display of images (SA, HLA, and VLA) that 
have been auto-realigned about the long axis of the heart and the defect maps and 
extents have been recalculated; and 

Figure 9 depicts a graph showing the mean change in the perfusion 
estimates when the auto-realignment is employed. 
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DETAILED DESCRIPTION OF EMBODIMENTS OF 
THE PRESENT INVENTION 

As required, detailed embodiments of the present invention are 
5 disclosed herein. However, it is to be understood that the disclosed embodiments 
are merely exemplary of an invention that may be embodied in various and 
alternative forms. Therefore, specific functional details disclosed herein are not to 
be interpreted as limiting, but merely as a representative basis for the claims and/or 
as a representative basis for teaching one skilled in the art to variously employ the 
10 present invention. Further, the image processing steps described below can be 
carried out automatically, using a computer. It is to be understood that persons 
skilled in the art understand how to program a computer to carry out these steps. 

Figure 3 depicts a flowchart of a reorienting process according to one 
embodiment of the present invention. As depicted in block 18 of Figure 3, ECT 
15 image data for a patient is acquired during a myocardial ECT study of the heart. In 
certain embodiments, patients used in the study may have known or suspected 
coronary heart disease. Separate ECT images are acquired with the patient at (i) rest 
and at (ii) stress (e.g. during exercise). 

In certain embodiments, the patient can undergo a stress Sestamibi 
20 of 0.6 GBq to rest Tetrofosmin of 1.3 GBq study for acquiring data with a 
simultaneous transmission/emission (TCT/ECT) system. Such a system is disclosed 
by Ficaro, et al., in "Simultaneous Transmission/Emission Myocardial Perfusion 
Tomography", Circulation, 1996; 93:463-473, which is incorporated herein by 
reference. The following parameters can be used with the simultaneous TCT/ECT 
25 system: 360 degree orbit, 60 steps per detector, 16 seconds per step, 64 by 64 
matrix size, and 6.3 millimeter pixel size. In certain embodiments, gated data can 
be reconstructed from 180 degrees projection image orbits from right anterior 
oblique (RAO) to left posterior oblique (LPO) using filtered back-projections. 
Transverse images can be filtered using a Butterworth filter with an order of about 
30 5.0 and a cut-off of about 0.30. 
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As depicted in block 20 of Figure 3, the initial center and axial limits 
of the left ventricle of the heart are obtained. Turning to Figure 2a, the axial limits 
are represented by slice positions 5 and 25 for the apex limit (z apex ) and base limit 
(z base ), respectively. While a user can input these values, the values can also be 
5 generated automatically from ECT images of the LV. For instance, to find x c and 
y c , a summed image from the axial and mid-third of the short axis (SA) image 
volume is made. The initial estimate for x c is the first moment of a linegram 
constructed from the upper half of the summed image. The initial estimate for y c 
is the first moment of a linegram constructed from the right 5/8 of this summed 

10 image (columns 24-64 for a 64x64 image matrix). Using these seeds as the center 
of an 8x8 neighborhood, the final estimates for x c and y c are assigned by fitting an 
ellipse to the radial information from the circumferential profile analysis. For the 
axial estimates (z c , z apex , z base ), a linegram involving only the anterior portion of the 
mid-axial slice down to the y c estimate is used. A threshold of 55% of maximum 

15 activity is applied to the linegram, from which the z apex and z^ are determined. The 
value for z c is the integral midpoint between z apex and z base . 

In block 22, the endocardial and epicardial surfaces of the left 
ventricle are determined. In certain embodiments, an LV surface detector can use 
a two-dimensional gradient image to determine initial estimates for the endocardial 

20 and epicardial surfaces. Weights based on activity, motion, and consistency of the 
endo- and epicardial surface estimates in space and time are assigned. A series of 
one- and two-dimensional weighted splines are used to refine the surfaces estimates. 
Since splines insure the continuity of the first derivative, the natural curvature of the 
myocardium is preserved even in the presence of large and severe perfusion defects. 

25 Using the new surface estimates, threshold values as a function of position are 
determined and used in conjunction with a morphological operator to create a 
segmented image of the LV volume. The segmented image is used for contiguity 
information in updating the spline weights, and the spline interpolators are applied 
again to better refine the surface estimates utilizing the contiguity information. With 

30 the new surface estimates, the LV wall activity is fitted to a Gaussian function to 
estimate the mid- wall position and thickness. For gated studies, the volume curve 
is estimated to determine the end-diastolic (ED) and end-systolic (ES) frames. The 
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average FWHM value for ED is used to determine the physical wall dimensions and 
is scaled to be an average thickness of 10mm. Wall thickness estimates for the 
remaining frames in the R-R cycle are adjusted based on the wall thickness scale 
factor used at ED, and myocardial mass is conserved. For ungated studies, the 
5 average wall thickness is assigned to be 15mm. The spline interpolating functions 
are used a final time to provide smooth estimates for the endocaridal, epicardial and 
wall thickness estimates. 

In block 24, the slice range for reorientation is determined. A slice 
is a two-dimensional SA image of the LV. The slice range is a portion of the LV 

10 along the total length (L) in which the LV is a shape similar to a donut. Turning to 
Figure 2a, the total length is represented by points 5 and 25, representing L=21 
slices. The axial center of the LV is defined as slice position 15. In certain 
embodiments, the slice range can be about 45% of the LV length centered at the 
midpoint of the LV. In other words, the slice range is essentially 0. 15 times L from 

15 the base to 0.60 times L from the base. For the study depicted in Figures 2a-2c, the 
slice range is from slices 12 to 22. This number of slices, 11, is defined as N. This 
slice range is typically selected because the LV is shaped similar to a cylinder and 
the endocardial and epicardial surfaces are elliptical in shape. 

In block 26, center coordinates are determined for each slice i in N 
20 as defined in block 24. Using the endocardial surfaces, i.e. inner surfaces of the 
heart, the center coordinates x[i= 1 to N], y[i= 1 to N] for each SA slice in i = 1 to 
N are determined from a fit of the surface points to an ellipse. In certain 
embodiments, the surface points are greater than or equal to 60. 

According to block 28, translation (A) values and rotation (0) values 
25 are determined. From the array of center coordinates, the rotation values 0 X , 0 y and 
translation values Ax, Ay along the x and y axes required to reorient the long axis 
of the heart to the z-axis can be determined. In certain embodiments, extreme 
values, defined as 6 X , 0 y greater than about 10 degrees for rotation values or Ax, Ay 
greater than 10 pixels for translation values, are dampened during this step of the 
30 process. 
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According to decision block 30, the process decides whether 
reorientation is necessary. In certain embodiments, any translation (A) value greater 
than about 0.1 pixels or any rotation (0) value of greater than 0.1 degrees 
necessitates reorientation. If reorientation is necessary, the LV is reoriented, as 
5 depicted in block 32. The reorientation can include multiple steps. In certain 
embodiments, the first step is determining a reorientation matrix for reorienting the 
current LV volume based on 0 X , 0 y and Ax, Ay. In certain embodiments, the second 
step is updating the cumulative reorientation matrix for reorienting the original LV 
volume. In certain embodiments, the third step includes updating the LV center by 

10 transforming original LV center position by using the cumulative reorientation 
matrix. These transformation steps can be performed by using an affine transform, 
although other transforms known in the art can also be utilized. Once the 
transformation is carried out, the steps described in blocks 22, 24, 26, and 28 are 
carried out until reorientation is not necessary. Once reorientation is not necessary, 

15 the cumulative reorientation matrix is stored, as depicted in block 34. According 
to block 36, functional and perfusion parameters can be determined for the 
reoriented ECT images. Following is a non-limiting example of the beneficial 
results that can be obtained using the process identified in Figure 3. 

EXAMPLE 

20 The patient population included thirty -one (31) consecutive patients 

with known or suspected coronary heart disease with no exclusions based on body 
habitus. Each patient underwent a stress Sestamibi of 0.6 GBq versus rest 
Tetrofosmin of 1.3 GBq study acquired with a simultaneous TCT/ECT system. The 
acquisition protocol included, but was not limited to, the following parameters: 360 

25 degrees orbit, 60 steps per detector, 16 seconds per step, 64 x 64 matrix, 6.3 mm 
pixel size. Further, all gated data were reconstructed from 180 degrees projection 
image orbits from RAO to LPO using filtered back projection. Transverse images 
were filtered using a Butterworth filter (order 5.0, cut-off 0.30). 

For each of the stress gated studies, three short axis (SA) volumes 
30 were constructed. Visually ideal alignment to the long axis of the heart, misaligned 
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+/- 5 degrees in the transaxial and sagittal oblique planes, misaligned +/- 10 
degrees in the transaxial and sagittal oblique planes were constructed. A single SA 
volume was constructed for the rest study. 

The baseline for this example was no realignment. Functional and 
5 perfusion data were quantified using a software program, 4D-MSPECT, available 
from the University of Michigan (Ann Arbor, Michigan) which did not adjust for 
misalignment of the SA volume about the long axis. 

Realigned image volumes were constructed using an automatic 
computer-implemented process of the present invention. Functional and perfusion 
10 data were quantified using the realigned image volumes. 

The following parameters were investigated: left ventricular ejection 
fraction (EF), summed stress score (SSS) for perfusion, total perfusion defect extent, 
and total perfusion ischemic extent. For the misaligned (MA) and the auto-aligned 
(A A) (according to the present invention) data sets, maximal differences between the 

15 3 stress data sets were recorded. Significant differences between the MA and A A 
data sets were determined using the Students t-test. Values of p less than about 0.05 
were considered significant. Values of p less than about 0.05 were considered 
significant. The t-test determines the degree of correlation between two discrete 
variables. When the p-value, the result of this test, is less than 0.05, the probability 

20 is greater than a 95% that the MA and AA distributions are different. 

Figure 4 depicts a splash display of images (SA, HLA, and VLA) 
with various rotational misalignments. The variability in EF values between the 
studies is substantial, that is, 41%, 39%, and 47%. 

25 Figure 5 depicts a splash display of images (SA, HLA, and VLA) that 

have been auto-realigned about the long axis of the heart. The variability in EF 
values was substantially decreased relative to the EF values of Figure 4, that is, 
45%, 45%, and 43%. 
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Figure 6 depicts a graph showing the mean change in the functional 
estimates when the auto-realignment is employed whereas the mean maximal 
difference percentage for EF based on the MA images was 2.54 _+ 1.62, the value 
for the AA images was 1.01 +_ 0.72. 

Figure 7 shows a splash display of images (SA, HLA, and VLA) with 
various rotational misalignments and the corresponding defect blackout maps and 
regional defect extent values. For example, the SA G-Str 5 defect blackout may 
include vascular regions with 7% blackout, 1% blackout and 1% blackout. The 
blackout regions may cause a false diagnosis for coronary heart disease. 

Figure 8 shows a splash of images (SA, HLA, and VLA) that have 
been auto-realigned about the long axis of the heart and the defect maps and extents 
have been recalculated. As can be seen by Figure 8, the blackout maps show 0% 
blackout. The AA images result in less false diagnoses relative to the MA images 
of Figure 7. 

Figure 9 depicts a graph showing the mean change in the perfusion 
estimates when the auto-realignment is employed. The y-axis is the mean maximal 
differences percentage and the x axis is variability in perfusion estimated using MA 
and AA images for SSS, defect extent, and ischemia extent. 

In addition to the results graphed in Figures 6 and 9, the 
automatically realigned image sets displayed no detectable differences compared to 
visual ideally reoriented images irrespective of the magnitude (less than about 10 
degrees) and angles of misalignment. Further, regional differences by individual 
vascular territory were even greater and of clinically significant magnitudes 
compared to the global measures reported in Figures 6 and 9. Advantageously, 
means and standard deviations for each of the compared parameters were 
significantly reduced by automatic alignment including centering. 

While embodiments of the invention have been illustrated and 
described, it is not intended that these embodiments illustrate and describe all 
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possible forms of the invention. Rather, the words used in the specification are 
words of description rather than limitation, and it is understood that various changes 
may be made without departing from the spirit and scope of the invention. 
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